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ABSTRACT 

The X-ray source IGR J16318— 4848 was the first source discovered by IN- 
TEGRAL on 2003, January 29. The high energy spectrum exhibits such a high 
column density that the source is undetectable in X-rays below 2 keV. On 2003, 
February, 23-25, we triggered our Target of Opportunity (ToO) program using 
the EMMI and SOFI instruments on the New Technology Telescope of the Euro- 
pean Southern Observatory (La Silla) to obtain optical and near-infrared (NIR) 
observations. We confirm the already proposed NIR counterpart and for the 
first time extended detection into optical. We report here photometric measure- 
ments in the i?, J, J, bands, upper flux limits in the Bh and V bands, lower 
flux limits in the H and Kg bands. We also obtain NIR spectroscopy between 
0.95 and 2.52 /im, revealing a large number of emission lines, including forbidden 
iron lines and P-Cygni profiles, and showing a strong similarity with CI Cam, 
another strongly absorbed source. Together with the Spectral Energy Distribu- 
tion (SED), these data point to a high luminosity, high temperature source, with 
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an intrinsic optical-NIR absorption greater than the interstellar absorption, but 
two orders of magnitude below the X-ray absorption. We propose the following 
picture to match the data : the source is a High Mass X-ray binary (HMXB) at 
a distance between 0.9 and 6.2 kpc, the optical/NIR counterpart corresponds to 
the mass donor, which is an early-type star, maybe a sgB[e] star, surrounded by a 
dense and absorbing circumstellar material. This would make the second HMXB 
with a sgB[e] star as the mass donor after CI Cam. Such sources may represent 
a different evolutionary state of X-ray binaries previously undetected with the 
lower energy space telescopes ; if it is so, a new class of strongly absorbed X-ray 
binaries is being unveiled by INTEGRAL. 

Subject headings: stars: circumstellar matter — stars: emission-line. Be — X-rays: 
binaries 



1. INTRODUCTION 

ICR J16318-4848 was the first new source discovered by the INTEGRAL IBIS/ISGRI 
imager (Ubertini et al. 2003; Lebrun et al. 2003). The source was first detected on 2003 
January 29 (Courvoisier et al. 2003) in the energy band 15-40 keV with a mean 20-50 keV 
flux of 6 X 10^^ ergcm~^s~^. The source was thereafter regularly observed for two months. 
No X-ray counterpart could be found in the ROSAT All Sky Survey (Voges et al. 1999). The 
source was observed by XMM-Newton on 2003 February 10, which detected a single X-ray 
source within the INTEGRAL error box using the EPIC PN and MOS cameras (Striider et 
al. 2001; Turner ct al. 2001), giving the most precise position to date : a = 16'^31™48';6, 

5 = -48°49'00" with a 4" error box (a circle with a radius of 4")(Schartcl ct al. 2003). 
X-ray spectroscopy revealed a very high absorption column density : iVn > 10^"^ cm^^ (Matt 

6 Guainazzi 2003; Walter et al. 2003), which renders the source invisible below 2 keV. This 
amount of absorption is unusual in Galactic sources. This could explain the non detection 
by ROSAT, although the source was discovered at a similar flux level in archival ASCA 
observations in 1994 (Murakami et al. 2003; Revnivtsev et al. 2003) on both CIS and 
SIS instruments (between 0.4 and 10 keV). Relatively bright and highly absorbed sources 
hke ICR J16318— 4848 could have escaped detection in past X-ray surveys and could still 
contribute signiflcantly to the Galactic hard X-ray background in the 10-200 keV band. 



^ Based on observations collected at the European Southern Observatory, Chile (observing proposal ESO 
N 70.D-0340). 
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The high column density prompted a counterpart research in the near-infrared (NIR) : in 
the EPIC error box, a possible counterpart was proposed by Foschini et al. (2003) using the 
Two Micron All Sky Survey (2MASS) with the following magnitudes : J = 10.162 ± 0.018, 
H = 8.328 ± 0.037, K,, = 7.187 ± 0.015 (Walter ct al. 2003). On the other hand, no radio 
emission at the position of the source could be detected. These broad characteristics suggest 
that it is an X-ray binary, the nature of the compact object (neutron star or black hole) 
being subject to debate. 

In the course of a Target of Opportunity (ToO) program at the European Southern Obser- 
vatory (ESO) dedicated to look for counterparts of high energy sources newly discovered by 
satellites including INTEGRAL (PI S. Chaty), we carried out photometric observations in 
the optical and NIR, and spectroscopic observations in the NIR. The goals were to search for 
likely counterparts within the EPIC error box, to obtain informations about the environment 
and the nature of the source, especially about the mass donor. 

The main results of this paper have been already presented in a very condensed way in Chaty 
& Filliatre (2004). In section 2, we describe briefly our observations and their reduction ; in 
section 3, we report on astrometry and photometry of all possible counterparts of the source, 
and give the most likely candidate ; for the latter, we give in section 4 the spectral energy 
distribution, putting our data together with the survey archives and published data, and 
derive absorption and temperature estimates ; in section 5, we study the spectral lines ; in 
section 6, we discuss the distance to the source, its nature and the nature of its components ; 
then, in section 7, we conclude. 

2. THE OBSERVATIONS 

The observations were carried out using the New Technology Telescope (NTT) of the 
ESO at La Silla (Chile), on 2003 February 23-25. We used the EMMI (optical) and SOFI 
(NIR) instruments. For astrometry and photometry, we had a total of 8 broad bands : Bb, 
V, R, I, Z, J, H, Ks- For spectroscopy, we used 2 bands : GBF (Grism Blue Filter, range 
0.95 — 1.64 iim, resolution 930) and GRF (Grism Red Filter, range 1.53 — 2.52 /^m, resolution 
980). 

All the data were reduced using the IRAF suite^. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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2.1. Optical Photometric Data 

The optical data have been taken on 2003 February 24 with EMMI between UT 7h30 
and 8h30. The seeing was typically 0'.'9, and the airmass on the source was between 1.225 and 
1.277. We acquired one 180 second exposure for each band {Bb, V, R, I, Z) and an extra 60 
second exposure for Bb. The data analysis involved subtraction of the bias and division by the 
flat field. For calibration, we obtained an exposure of 30, 15, 10, 10, 10 seconds, respectively 
of the five stars of PG1633+099 of Landolt catalog (Landolt 1992). The EMMI detector 
is a mosaic of two arrays, each with two readouts, resulting in a 4152 x 4110 array with a 
pixel-scale of 0'.'166 per pixel. This makes four parts with slightly different characteristics. 
The source fell into the second readout of the first array, as did three of the standard stars, 
therefore we used these three stars to compute the color equation for BbVRI. These stars 
were PG 1633+099 B, C and D. As the airmass was identical for these three stars, it was 
then not possible to determine the extinction coefficient with them, and we adopted the 
values reported in the EMMI manual (Gonzalez et al. 2002). We found the zero-point and 
the color coefficient by a hnear regression. The results are given in Table 1. The zero-point 
magnitude is in good agreement with data given in Gonzalez et al. (2002) (converted for 
1 ADU per second to be directly comparable with our results) : 24.98 ± 0.03, 25.69 ± 0.01, 
25.98 ± 0.02 and 25.28 ± 0.03 for BbVRI, respectively. 



2.2. NIR Photometric Data 

The NIR photometric data have been taken with SOFI on the NTT between UT 7 and 
10 h on 2003 February 23. The SOFI detector is a 1024 x 1024 CCD with a pixel-scale of 
0'.'288 per pixel. The seeing was typically 1", and the airmass on the source was between 
1.088 and 1.388. We obtained 9 exposures of 10 seconds and 11 exposures of 5 seconds in 
the J band. We obtained 5 exposures of 10 seconds and 9 exposures of 2 seconds in the 
H band. We acquired 245 exposures of 2 seconds in the Ks band. In order to correct for 
the bright NIR sky, following the standard procedure, the pointings of all these frames were 
dithered. This allowed us to obtain a sky template by median filtering the frames. This 
template was then subtracted from the flat-fielded frames. The ahgnment of the frames was 
performed using ~ 2000 common stars. Within each band, frames with different exposure 
times were averaged separately. For calibration, five exposures for each band of standard 
stars sj9157 (integration time 2 seconds) and sj9170 (1.2 seconds) of Persson et al. (1998) 
were done. The color term was taken from Lidman et al. (2002), because we had three 
bands and only two standard stars. We found the zero-point and the extinction coefficient 
by a linear regression. The results are reported in Table 1. The zero-point magnitude is in 
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good agreement with data given in Lidman et al. (2002) : 23.2, 23.1 and 22.5 ± 0.01 for J, 
H and Kg, respectively. 

2.3. NIR Spectroscopy 

The NIR spectroscopic data have been taken with SOFT low resolution grism on the 
NTT between UT 8 and 9 h on 2003, February, 25 with the GBF and GRF filters. The 
seeing was typically 0'/9, and the airmass on the object varied between 1.160 and 1.260. For 
each GBF and GRF bands, wc obtained 24 exposures of 6 seconds. After subtraction of the 
bias, the exposures were divided by the response of the fiat field along the dispersion line. 
The resulting individual rough spectra show a dispersion of around 20 %, which seems to 
be random, and not associated to a short time variation of the source. Therefore, the two 
averaged spectra have a relative precision of 4 % in flux. Given the seeing and the slit width 
of 1", we expect to collect around 75% of the flux ; however, as we will scale the spectra 
to match the broad band magnitudes, absolute calibration is not an issue. Wavelength 
cahbration was done by observing a Xenon arc. We obtained 16 exposures of 5 seconds 
for standard star hip80456 (my = 7.74, spectral type F5V), 8 exposures of 10 seconds 
for standard star hip83612 (my = 8.38, spectral type GIV). With these stars we made a 
telluric correction, followed by a flux calibration using a blackbody spectrum extrapolated 
from the magnitude and spectral type of the standard star. Some absorption lines especially 
Br7, were seen on the standard stars : we removed the lines and replaced them by a linear 
interpolation, as the continuum was smooth enough. As for the object the two bands did 
not match in continuum flux in the overlap region, the continuum was extracted and scaled 
to match the J and Kg broadband values, giving a good agreement with the H broadband 
value (see subsection 3.2). In two regions, between 1.35 and 1.45 /im, and between 1.8 and 
2.0 /im, there is a strong absorption and the flux calibration did not work successfully. The 
spectrum is also very noisy above 2.35 /im. 

3. IDENTIFICATION OF THE COUNTERPART 

3.1. Astrometry of the candidates 

The EMMl pixel-scale is 0'.'166. Therefore, the error box of 4" reported by EPIC (Schar- 
tel et al. 2003) corresponds to a circle with a radius of 24 pixels, noticeably greater than 
the FWHM~ 5 pixels measured on the stars on the processed BbVRIZ frames, and corre- 
sponding to the 0'.'9 seeing. To find the pixel coordinates of the center of the EPIC error 
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box, we proceed as follows : 

• we use the GSC^ catalog to find the astromctric coordinates of 15 stars in the field 
with a precision compatible with the EMMI pixel scale ; 

• we use the relevant frame of the DSS catalog^ to find visually these stars and mark 
them on the EMMI frames ; this DSS frame is not actually used for astrometry, only 
for visual identification ; 

• we compute the transformation matrix to go from the astrometric coordinates to EMMI 
pixel coordinates ; the residuals are well below the error box, at the 0'/2 level, compatible 
with the EMMI pixel-scale and the 0f!5 error reported for GSC^ ; 

• we use this matrix to compute the pixel coordinates of the center of the error box. 

We show the position of the source with the EPIC error box of 4" superimposed, on EMMI 
and SOFI frames, respectively on Figure 1 and Figure 2. 

Given the SOFI pixel-scale, the uncertainty of 4" reported by EPIC corresponds to a circle 
with a radius of 13.9 pixels, noticeably greater than the fwhm~ 3 pixels measured on the 
stars on the processed JHKg frames. We use a similar method to find the pixel coordinates 
of the center of the EPIC error box, but we use 2MASS^ as catalog. Again, the residuals 
in the construction of the transformation matrix are about a pixel size, e. g. (T.'S, with an 
uncertainty for the 2MASS catalog about (y/OQ. 

On Figures 1 and 2, the scale is nearly the same, and given by the size of the 4" error box, to 
allow for visual comparisons. The astrometry is consistent within the errors in both optical 
and NIR image sets. Three candidates appear within the error box. For further references, 
we will label them 1, 2 and 3 on Figure 2, left panel. The position of the center of each 
candidate with respect to the center of the circle is, on frame J : 3" SW. 3'.'6 NE and 2'.'2 
NE, respectively. Figure 2 also labels two field stars (4 and 5) for comparisons. 
We give on Table 2 the astrometry of the three candidates, computed by inverting the 
transformation matrix from GSC coordinates to pixels on frames /, J, and averaging, 
putting statistical weight according to the uncertainties on the catalogs ; the dispersion being 
below O'.'l, and given the uncertainties on GSC and 2MASS catalogs, we take the residual 
value of (y.'5 as error estimate. 

^ http : / /archive .eso.org/gsc/gsc 
^http: / /archive. eso.org/ dss/dss 

^http: / /www-gsss. stsci.edu/ gsc/ gsc2/GSC2home.htm 
^ http : / /ir sa . ipac . caltech . edu/ 
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Table 1. Terms of color equation. The extinction coefficients for the optical bands are 
taken from Gonzalez et al. (2002), the color terms for the NIR are taken from Lidman et 
al. (2002), the other coefficients are measured values. 



Band 


Zero-point mag. 


color term 


Extinction coeff. 


Bb 


25.03 ± 0.03 


-0.03 ± 0.01(5 -y) 


0.214^ 


V 


25.59 ±0.02 


0.08±0.01(K-i?) 


0.125^ 


R 


25.90 ±0.03 


-0.005 ±0.2(7? -J) 


0.09P 


I 


25.42 ±0.04 


0.04±0.29(i?-/) 


0.05P 


J 


23.34 ±0.06 


-Q.m{J -Kf 


0.40 ±0.05 


H 


23.17 ±0.06 


-Qm2{J - Kf 


0.31 ±0.05 


Ks 


22.39 ± 0.06 


0.023(J - Kf 


0.15 ±0.05 



^Not measured, taken from Gonzalez et al. (2002) 
•^Not measured, taken from Lidman et al. (2002) 




is left. The scale is given by the error box. Candidate 1 is at the south-west border of the 
circle, and candidates 2 and 3 are at the north-east border (see Figure 2 for identification). 
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Fig. 2. — The error box of 4" is shown on the J (shortest exposures), H (shortest exposures) 
and Ks frames obtained with SOFI. North is up, east is left. Star labels are given on the J 
frame. 




Fig. 3.— The USNO Bl.O plate, IIIaF emulsion, filter RG630. The EPIC error box is 
reported. North is up, east is left. 



Table 2. Astrometry of the three candidates, with an uncertainty of ±(y.'2. 



Candidate 


a 


5 


1 


16^31'"48'!31 


-48°49'00'.'7 


2 


16'^31°^48!89 


-48°48'57'^7 


3 


16''31™48!78 


-48°48'58'/7 
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3.2. Broad Band magnitudes of the candidates 

We performed photometry in the BhV RIJHKa bands using the daophot package of 
IRAF, for the three candidates reported in Figure 2. For Bb and V frames, we indicate for 
candidates 1 and 3 the magnitude of the faintest star as an upper hmit for the flux. The 
results are given in Table 3. Table 4 is a summary of archival observations of candidate 1. 
Candidate 2 is barely detectable in Bb, more clearly in V. Candidates 1 and 3 are detected 
in the other bands. Candidate 1 is also found in the 2MASS catalog and was suggested as 
an infrared counterpart by Foschini et al. (2003) ; it is saturating the EMMI detector in Z. 
It is also one of the brightest stars in the NIR field. It saturates in H and Ks (even with the 
shortest exposure time), and in J (shortest exposures) it peaks to 32000 ADU (Analogue 
to Digital Units), above the limit of linearity of 10000 ADU of SOFI (Lidman et al. 2002). 
The position is taken near the maximum, before saturation. Therefore magnitudes for this 
star reported in Table 3 are upper limits in H and Kg, and there is a possible bias in J ; 
note also that the optical frames show that there is no neighbour in the glare of the star, so 
it is likely that there is no significant contamination in NIR. 

Candidate 1 has been seen in the DENIS catalog^, with magnitudes in agreement with our 
results. It has been seen in the DSS2^ infrared plate (1980), but not seen in the DSS2 red 
plate (1992). However, as stated by Foschini et al. (2003), this candidate is not visible on 
the USNO Bl.O plates^ taken with IIIaF emulsion, filter RG630 (see Figure 3) and OG 590, 
corresponding to the R band, obtained respectively in 1982 and 1992. This is in contradiction 
with Walter et al. (2003), who identified a counterpart in the USNO Bl.O catalog ; however, 
it appears at the light of our data that their identification corresponds to our candidate 
2, whereas their identification with 2MASS is actually our candidate 1. This mistake was 
probably due to the fact that candidates 1 and 2 have nearly the same magnitude in R in 
our results, and that candidate 2 is reported to vary of one magnitude over an interval of 50 
years (between the two epochs of the USNO catalog). Our observation of candidate 1 in the 
R band is therefore, up to our knowledge, a discovery. As the magnitude limit on the USNO 
Bl.O catalog is around 20, our detection seems to imply a variability in R band of more 
than two magnitudes, whereas there is no variability on J between 2003 (our data) and 1999 
(2MASS). Such a strong variability in the optical (although comparisons with photographic 
measurements of objects of such unusual colors may be rather uncertain, because of slightly 
different bandpasses), associated to the high variability of ICR J16318— 4848 in the X-rays, 



''http://cdsweb.u-strasbg.fr/denis.html 

^http: / /archive. eso.org/ dss/dss 

^http: / / www.nofs.navy.mil/data/fchpix/ 
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would suggest that candidate 1 is the genuine counterpart of IGR J16318— 4848. 

We also seek for the source in the IRAS Point source catalog V2.1^°, and find no source 
closer than 230". We found 1900 IRAS sources in a 3° radius around the EPIC position. We 
choose as a robust upper limit for Table 4 the flux such as 90% of the IRAS sources have 
greater fluxes : indeed, the flux distribution shows a obvious deflcit of fainter sources. 
Radio observations were performed with the Australia Telescope Compact Array (ATCA) 
at 4.8 and 8.6 GHz (Walter et al. 2003) on 2003, February 9, and show no detection with a 
1 a upper limit of 0.1 mJy. 

The differences in magnitude R — J, R — Kg, J — H, J — Ks of the stars of the field are 
shown on Figure 4. The unusual colors of candidate 1 are obvious, suggesting a reddening 
by absorption : the star is within the bulk of the distribution in the J — H, J — Kg plane, 
but is a clear outlier in the R — J, R — Kg plane, where absorption is more critical. As being 
brighter at longer wavelength is expected for a highly obscured source, this suggests that 
candidates 2 and 3 may be field stars, and candidate 1 is therefore the most likely candidate 
for IGR J16318-4848 counterpart. 



3.3. Summciry 

In the two previous subsections, we have shown that three sources have positions com- 
patible with the error box given by EPIC. Candidate 1 is seen on 2MASS, and we confirm 
its photometry in the J band. We discovered this candidate in the R band, since it is 
not detected in the USNO-Bl.O catalog. The dependence of its magnitudes on the wave- 
length shows a strong brightening at longer wavelengths, which indicates that it is strongly 
absorbed, whereas the two other candidates have behaviors comparable to the field stars. 
Therefore, in the rest of this paper, we will consider candidate 1 as the genuine counterpart 
of IGR J16318-4848. 



4. THE SPECTRAL ENERGY DISTRIBUTION 

4.1. Spectral continuum 

The continua of both GBF and GRF spectra were extracted, and scaled to match our 
J magnitude and the 2MASS Kg magnitude for candidate 1. This method is approximate. 



'http://irsa.ipac.caltech.edu/IRASdocs/iras.html 
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Table 3. Photometry of the stars labeled in Figure 2. 



Band 


1 


2 


3 


Bb 


> 25.4 ± 1 


22.70 ±0.12 


> 25.4 ± 1 


V 


> 21.1 ±0.1 


19.67 ±0.03 


> 21.1 ±0.1 


R 


17.72 ±0.12 


17.86 ± 0.05 


20.25 ±0.14 


I 


16.05 ±0.51 


17.71 ± 0.06 


19.71 ±0.17 


J 


10.33 ±0.14 


16.42 ±0.14 


16.97 ±0.14 


H 


< 10.35 ±0.15 


16.43 ±0.16 


16.93 ±0.16 


Ks 


< 9.13 ±0.10 


15.22 ±0.20 


14.86 ±0.17 



Table 4. Archival data on candidate 1. 



Name 


Year 


Band 


Result 


DSS2 


1980 


Infrared plate 


seen 




1992 


Red plate 


unseen 


USNO Bl.O 


1982 


RG630 


unseen 




1992 


OG590 


unseen 


2MASS 


1999 


J 


10.162 ±0.018 






H 


8.328 ±0.037 






Ks 


7.187 ±0.015 


DENIS 


2003 


I 


16.217 ±0.06 






J 


10.239 ±0.05 






K 


7.255 ±0.07 


IRAS 


1983 


12 /xm 


< 0.9 Jy 






25 //m 


< 0.7 Jy 






60 /.tm 


<4Jy 






100 //m 


< 40 Jy 


ATCA 


2003 


8.6 GHz 


< 0.1 mJy 






4.8 GHz 


< 0.1 mJy 
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Fig. 4. — Color of stars of the field, detected together in R, J, H, and Kg. 2MASS magnitudes 
were used in H and Kg for candidate 1. The numbers 1 to 5 refer to stars lebelled on Figure 2. 
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given the unusual colors of the object, nevertheless this leads to a very good agreement in 
the overlapping region, and with the corresponding 2MASS H magnitude. This approach 
assumes that the observed magnitudes are mainly due to the contimia ; indeed, both in 
GRF and GBF, the power of the continua represents about 95 % of the total power. In 
Figure 5, we show the Spectral Energy Distribution (SED) in (z/, vF{v)) units, including 
the two continua, with results from Table 3 and the archival data of Table 4. At this 
stage, no dereddening is done. We have also included in the SED a compilation of several 
observations of the source, covering a very wide range of wavelengths, going from centimeter 
(ATCA) to X-rays (XMM, INTEGRAL), over 10 decades in wavelength. Data confirming 
our observations, as for instance the DENIS catalog, are not reported on the SED. However, 
two caveats have to be mentioned : 



• the data were taken at different dates, and the source is variable (sec Walter et 
al. (2003), and our comment on the non-detection of the source on the USNO Bl.O 
catalog) ; 

• the emissions at different wavelengths can come from different components of the sys- 
tem, as the source is likely to be a complex system involving a compact object, a 
companion star of an early type, and dust (see below). 

Therefore, throughout the discussion, we will emphasize on our optical/NIR data, and use 
the other data mainly for completion or confirmation. 



4.2. Absorption and temperature estimates 

We suppose that our NIR (photometry and spectroscopy) and optical data are mea- 
surements of a blackbody at temperature T with absorption Ay. If the emitting region is 
spherical of radius r at a distance D from Earth, then the measured flux in mJy is given by : 

where A^, is the absorption at frequency z/, given Ay, and is computed according to the 
formulae of Cardelli et al. (1989). We fit our R, 7, GBF and GRF fluxes for parameters 
r/D, T and Ay using a gradient minimization method. The fit shown on Figure 5 has 
r/D = b ■ 10~^°, T = 20 250 K, Ay — 17.5, but because of a strong degeneracy, this is not 
the sole choice for these parameters. 

This fit reproduces our data well, rendering plausible the hypothesis that an absorbed black 
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Fig. 5. — Observed SED in [u, vF{y)) units, including the results of Table 3 and the archival 
data of Table 4. At this stage, the rescaled continuum spectra in GBF and GRF and 
literature data. The B and V data are upper limit only. The dashed curve corresponds to 
the fit discussed in section 4.2 : an absorbed black body, representing well the data. The 
results of INTEGRAL, XMM, IRAS and ATCA are also shown. 
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body is the origin of the continuum emission. Note that the absorption law is unknown 

for wavelengths above 3.3 /xm (9 x 10^^ Hz) (Cardelli et al. 1989), and we assume that the 
expression is still valid for longer wavelengths. The value of the fit is quickly decreasing at 
frequencies lower than the Kg band, rendering likely the non detection of the source by IRAS 
and ATCA : at 2.5 x 10^^ Hz, the average flux of the fit is only 16 % of the upper limit of 
IRAS. 

However, the parameters of the fit are strongly degenerate, and therefore the parameter esti- 
mates returned by the minimization method depend strongly on the first guess. Physically, it 
means that a cool, close source with low absorption can fit the data as well as a hot, distant 
and absorbed one. By analogy to chaotic mechanic systems, we call "attractor" the locus in 
the parameter space of the various minimum wc can get by varying the first guess. To 
trace this attractor, we use the minimization method 5000 times, with uniformly distributed 
parameters over r/D = IQ-^^ - 10~^ T = 1 000 - 61 000 K, Ay = - 25. The result is 
shown on Figure 6. There is a clear attractor solution, all individual points corresponding to 
undistinguishable fits with nearly identical x^. In the temperature/absorption plane, almost 
all points converged towards a line well fitted by : 

4^ = 1 - - (2) 

with Avo = 18.68, Tq = 1200K. Note that there is a slight discrepancy for high tempera- 
tures, indicating that this simple empirical model is not fully satisfying. The To parameter 
corresponds to the black body temperature that fit our data if there is no absorption. The 
Avo parameter is the absorption limit when T —>■ oo, i.e. the maximum absorption compat- 
ible with our data, and the convergence towards this limit is rather quick. Two regions can 
be loosely distinguished, corresponding to two different physical situations : 



a low temperature (below ~ 6 000 K) where the absorption is very badly constrained 
between Ay — 5 and 15. This temperature is compatible with a main sequence dwarf 
star photosphere, a cool red giant, or even with hot dust if the temperature is not 
too much above 1000 K. On the other hand, the heating of dust can be non thermal 
(free- free for instance). This case corresponds also to a higher r/D ratio, indicating a 
rather close and/or large source. 

a high temperature (above ~ 10 000 K), where the absorption is rather strongly con- 
strained between Av = 16 and 18.7. This temperature is compatible with an early 
type star photosphere. However, because the NIR corresponds to the Raylcigh- Jeans 
region of the spectrum, the temperature has almost no effect on the shape of the 
spectrum, and is therefore very weakly constrained. 
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This ambiguity of fitting with an absorbed blackbody had been already noted by Revnivtsev 
et al. (2003), using archival observations. We essentially agree with their results. 
The absorption can be intrinsic, or along the line of sight (interstellar absorption), or both. 
In the case of interstellar absorption, an upper limit (valid for a source not too close) can be 
obtained with the A'^h tool of HEASARC^^, using the results of Dickey & Lockman (1990). 
Within 1° around the position of the source, seven data points are available, well distributed 
and well centered. By averaging, and putting weight in order to favor closer points of the 
source, we obtain an estimate of the H I column density of 2.11 ± 0.3 lO^^cm"^. Using the 
relation of Predehl & Schmitt (1995), we get an estimate for the interstellar absorption 
of Ay — 11.8 ± 1.6. If the temperature is above ~ 10 000 K (the "high temperature case" 
above), this seems to point to an absorption more than 3 a above interstellar absorption, 
and hence is an indication of a strong intrinisic absorption, and that, if the star responsible 
for optical/NIR radiation is this hot, then we see it through a rather dense circumstellar 
material. This indication is in agreement with the fact that the source, if it is the candidate 
1 reported on Figure 4, exhibits unusual colors, whereas the neighbour stars labeled 2 to 5 
have color not different to those of field stars. However, given the proximity of the source 
to the Galactic plane, the interstellar medium may be patchy, hence the derived value of 
interstellar absorption must be taken only as an estimate. 

A more detailed modclization of the SED should include the contribution of the dust, taking 
into account the distribution of the material, and its physical properties. However, valuable 
contraints would come from data in the mid-infrared range, where the dust contribution is 
expected to be noticeable. In the Kg band, the contribution of the dust is likely to be hard to 
distinguish, and indeed our fit is satisfactory in this band. That is the reason why we choose 
to keep a simple absorbed blackbody model, being not in contradiction with the possible 
presence of dust. 



5. THE SPECTROSCOPY 

In this section, we focus on the spectroscopic lines in order to distinguish between the 
"low temperature" and the "high temperature" case of Figure 6. We will also gather key 
information to assign a spectral type to the mass donor of the source. 



^ http://heasarc.gsfc. nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 
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Fig. 6. — The attractor solution for the minimization procedure to fit an absorbed blackbody 
to data plotted on Figure 5. Upper plot is for absorption, lower plot is for r/D. The points 
converge towards an attractor, that is fitted by a dash line on the upper plot. The remaining 
spreaded points indicate that the 1 % tolerance criterion required for convergence has not 
been reached in 50 iterations. See text for a discussion of the fit. 
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5.1. Line identification 

The full GBF and GRF spectra are presented on Figures 7, 8 and 9. The GRF spectrum 
is cut in two parts for a clearer view. The spectra display around 80 lines, most of them 
being emission lines. All lines selected here are detected above the 2cr level, and the relative 
uncertainty in flux is roughly 12 %. 

Line fitting is performed with the splot routine of IRAF, assuming a Gaussian shape. 
Tables 5, 6 and 7 (corresponding to Figures 7, 8 and 9, respectively) report the central 
wavelength for the identification and for the fit. peak flux, equivalent width and FWHM. A 
negative figure for width indicates that it is an emission feature, a positive figure indicates an 
absorption feature. To assess if the lines are real features of the object, we also extracted the 
spectrum on a different aperture (for GBF only). This spectrum is an order of magnitude 
below the object spectrum and clearly featureless (except in the strongly absorbed region 
between 1.35 and 1.45 yum), proving that the sky subtraction worked successfully in order to 
remove the sky lines. 

Identification is done using the UKIRT catalog^^, and the works of Morris et al. (1996), 
Clark et al. (1999), Bandyopadhyay et al. (1997). Only the strongest lines are fitted in the 
noisiest regions. A question mark indicates that the identification is not certain, or that the 
fitted value is very sensitive to continuum parameters given in input to splot, and therefore 
highly uncertain. Several identifications were often possible, given our resolution. Note that 
there is always an ambiguity between hydrogen and He II lines, i.e. a transition He II (2n- 
2m) has nearly identical wavelength as the (n-m) hydrogen transition. In any case, there 
are only two very weak He II transitions with odd levels (25-8 and 20-9), suggesting that the 
contribution of even levels to hydrogen lines is negligible. Therefore, this ambiguity is not 
indicated in the tables. 

In very few cases, there seems to be a doTible peak structure, the lower peak is then indicated 
by a "?-" in the tables (for instance, the Pa (6-3) at 1.0941 /im is preceded by a peak at 
1.0911 fim). The separation is around 400kms~^, similar to the value reported for Ci Cam 
(Clark et al. 1999), but the result is by far too dubious to lead to any firm conclusion. The 
FWHM values for the lines of Xenon arcs are 21 and 29 A, and we found that line broadening 
is not significant for the well-defined lines. The gap between the wavelength of the identified 
transitions and the fitted values is consistent within 1 a level with no systemic velocity of 
the object : cAA/A = -110 ± 130 km s^^ 

Four lines show a P-Cygni profile, indicating that there is a stellar wind of circumstellar 
material. All four are He 1 : 2s3S-2p3po (1.0832 /xm), 3p3p°-5d3D (1.1970 /xm), 2s^S-2p^po 
(2.0587 //m), 3p^P°— 4s^S (2. 1126 //m). The P-Cygni profile is evident in the two latter cases. 



littp://www.jach.hawaii.edu/JACpublic/UKIRT/astrononiy 
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more dubious in the former two. Putting more statistical weight on the two last lines, we 
infer a wind velocity of 410 ± 40kms~^, the uncertainty corresponding to the dispersion of 
the result using the four lines. No other line shows any clear P-Cygni profile. 
No significant time variation of the lines (in flux or in wavelength) was found, but we cannot 
give constraints on this because of the short total duration of the observations (10 minutes 
between the first and the last GBF spectrum). 

Many spectral hues are common to IGR J16318-4848 and CI Cam (Clark et al. 1999) : 
47 hues out of the 57 identified hues of IGR J16318— 4848 are seen in CI Cam between 1 
and 2.35 //m (Clark et al. 1999), and conversely 42 over 68 identified hues of CI Cam are 
seen in IGR J16318— 4848 ; hues after 2.35 /xm were not listed by Clark et al. (1999), and 
the CO bands are not seen in the source, but Pfund lines are common with CI Cam ; an 
unidentified feature at 1.9855 fim is common to both sources ; the 2-10 keV X-ray spectrum 
of IGR J16318— 4848 (Matt & Guainazzi 2003) is very similar to the X-ray spectrum of CI 
Cam (Boirin et al. 2002). These are hints that the two systems are physically similar. 

5.2. Hydrogen Lines 

The Paschen, Brackett and Pfund series are obvious on the spectra. They allow us to 
possibly estimate the absorption, using the Brackett decrement, independently of the SED 
fitting. The relative intensities of the hydrogen recombination lines have been calculated by 
Hummer & Storey (1987) for an electron temperature of 10^ K and an electron density of 
lO^cm"^ and considering Case B recombination. We compute the integrated fluxes of our 
Brackett and Paschen lines and divide by these relative intensities to obtain the absorption 
law. We fit separately the Brackett GBF lines and the Brackett GRF lines (all the lines) for 
the absorption Ay, and put a very low statistical weight on lines for which the intensity is 
uncertain, e.g. those at the edges of the GBF and GRF spectra, or within the noisiest regions. 
Results are given in Figure 10. GBF and GRF lines, considered separately, give somewhat 
different absorption values, 25.4 and 18.1 ; the joint fit gives a value of Ay = 18.2. However, 
in the fitting procedure, the minimum value of absorption is not very sharply defined. Note 
also that Br (7-4) (or Br7) has a smaller flux than expected with respect to the higher order 
Brackett lines, suggesting that the fines are optically thick or self-absorbed. Nevertheless, 
this method clearly favors the high temperature part of Figure 6, because it indicates that 
absorption is higher than 16. Therefore, it is likely that we are actually observing an early 
type hot star through a dense absorbing column, rather than a cool late star with lower 
absorption. 
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5.3. Helium Lines 

The upper state of He I 1.0832 /iin doublet, a strong line on our spectrum, can be popu- 
lated by recombination, or can be collisionally excited from the metastable 2s'^S state, hence 
its flux can be large even when the ionized fraction of helium is low (McGregor, Hyland, 
& Hillier 1988). Therefore, this line can come either from ionized or neutral regions. The 
2.0581 /im transition is among the strongest in our spectrum : this could be a hint of the 
presence of a high-density region close to the star, in this case the 584 A transition would be 
optically thick (McGregor, Hyland, & Hilher 1988), but becaiisc of the strong absorption, 
this transition will be very difficult to observe. Because two mechanisms can populate these 
transitions, the 1.0832 /xm and 2.0581 /xm lines are not reliable indicators of circumstellar 
conditions. 

We may estimate abundance of helium from the intensity ratios of He I 1.197yLim and 
1.7007 //m to Br (7-4). This is a lower limit for the helium abundance, because as the 
ionization energy of helium is greater than the one of hydrogen we have A'"(He"'")/N(H+) < 
N(He)/N(H), and because He I will not arise in the shielded Fe''" zones. According to Allen, 
Jones, & Hyland (1985), we have : 

^(He+) _ ^ /(l-7007^m) _ /(1.197/im) 
N{R+) ■ /(Br(7-4)) ' /(Br(7-4)) ^ ^ 

The result depends on the absorption Ay we choose to deredden the lines, as the line wave- 
lengths are remote from each other. The impact of absorption on the second ratio is higher, 
because of a greater distance between the two wavelengths, and because the line at 1.197/im 
has an absorbed flux much weaker than the one at 1.7007 /xm, there must exist an absorption 
value for which the two ratios are equal. This happens for Ay — 17.4 ± 2.5, the uncertainty 
coming from a relative uncertainty of 4 % on the flux lines. It is remarkable that this value 
is within the attractor values in the high temperature region of Figure 6 (the value of 17.4 
is associated with a temperature of 18 000 K, however the constraint on temperature is very 
loose, as it can be seen in Figure 6). The corresponding value is iV(He"'"/H"'") = 0.32 ± 0.04. 
The first ratio is less subject to absorption, and gives A^(He"'"/H"'") values ranging from 0.2 
to 0.5 for Av going from 10 to 25. In any case, this is a value far higher than the solar value 
of A'"(He)/N(H) = 0.067, indicating an evolved star. Note however that Br (7-4) is smaller 
than expected considering the higher order Brackett series, and that its intensity, according 
to Figure 10 may be underestimated up to a factor 10, causing an equivalent overestimation 
of the A^(He+/H+) ratio. 

We also consider the relative intensities of lines at 1.0832 /im (doublet), 1.1970 fim, 1.7007 //m 
and 2.0587 /um, and compare to predicted recombination intensities (see e.g. Allen, Jones, 
& Hyland (1985)). The results are given in Table 8. For dereddening, we take Ay — 17.4. 
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Note the concordance (better than 1 %) between the 1.1970 /^m and 1.7007 //m, another in- 
dication that dereddening is correct. We are in a case very similar to the one of rj Carinae, 
discussed in Allen, Jones, & Hyland (1985), with a line at 1.0832 /xm about 4 times weaker 
than expected, and the prominence of the line at 2.0587 /xm. Very interestingly, the ratios 
for CI Cam show a similar behavior (Clark et al. 1999), with a line at 1.0832 fim about 3 
times weaker than expected, and with a 2.0587 /^m line 7 times relatively more prominent 
than in our data. In the case of rj Carinae, these facts were explained by invoking a large 
optical depth for 1.0832 //m, and a significant optical thickness in the l^S — l^P 584.3 //m 
transition of He for 2.0587 /xm. 

Two He II transitions, at 1.6241 and 1.7717 /xm, are present, although very weak. The first 
one was also seen in CI Cam (Clark ct al. 1999). They indicate a high temperature able to 
ionize helium. Rough flux ratio with Fe II lines, using the solar abundances, leads to flux of 
photons able to ionize helium and iron in agreement with a temperature between 10^ and 
2 X 10^ K. 



5.4. Metallic Lines 

Several Fe II lines can be seen in our spectra. Among the strongest, 1.1126 /im and 
2.089 fim are also seen in early-type high luminosity stars (McGregor, Hyland, & Hillier 1988), 
and arise in highly excited levels which are probably populated by UV fluorescence. 
A peculiar emphasis must be put on the presence of six forbidden iron lines. Although they 
are rather weak, the 1.1567 /im line being for instance weaker than its neighbors of He I and 
N I, they provide a very important clue for the identification of the star. 
The doublet emission of Na I 4s^S — 4p^P° is seen, although not separated. The emission 
is most likely fluorescent, being pumped by 0.3303 /xm photons. The ionization potential is 
low (5.1 eV), so that sodium will be largely ionized around the star. Therefore, the sodium 
emission comes from regions that are not exposed directly to the radiation of the star or of 
the compact object. This suggests the presence of matter surrounding both objects. The 
other lines of Na I, such as 1.1385 and 1.1407/xm, are not seen. 

The Mg II doublet at 2.138 and 2.144 //m is clearly seen. It can be excited by L/? fluorescence. 
In our case, the 2.138 //m line is twice more intense than the 2.144 //m line, in agreement 
with the higher statistical weight of the 5P3/2 transition, and also to its shorter distance to 
L/3. 
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5.5. Moleculcir Lines 

No lines from H2 have been found, suggesting either that shock heating is relatively 
unimportant, or that temperature is sufficiently high to dissociate the molecules. Molecular 
hydrogen has not been seen either in CI Cam (Clark et al. 1999), nor in the massive stars 
in transition studied by Morris et al. (1996). Emission lines of H2 are more usual in massive 
young stellar objects (Kumar, Bachiller & Davis 2002), suggesting that there is no such 
radiation-protected medium in the environment of this source. 

No CO bands have been found, although a strong 8 bin rebinning of the GRF spectrum 
shows a step of lOmJy over 0.1 /xm at 2.295 ^um, locus of the (2-0) bandhead. However, 
this is too weak and too uncertain. The CO flux, if present, is therefore lower than ~ 10 mJy 
(with no dereddening) . 

5.6. Summary 

This study on the spectroscopic lines has shown that there is a dense circumstellar ma- 
terial around the source, possibly enshrouding both the mass donor and the compact object, 
and with a stellar wind component. This points to a high absorption, and we have seen 
that the value of Ay = 17.4 is favored from the data. The line ratios favor a high tempera- 
ture ; along with the high absorption, this is consistent with the case "high temperature" of 
Figure 6, and points to a hot early-type star observed through a dense absorbing material. 
Moreover, the presence of forbidden fines and the fact that almost all the lines are in emission 
point towards a "[e]" star, e.g. a B[e], with a striking similarity with the companion star in 
CI Cam. 




Fig. 7. — The GBF spectrum for the object. Indicated hnes refer to Table 5. 




Fig. 8. — The GRF spectrum for the object, first part. Indicated hnes refer to Table 6. 
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Fig. 9. — The GRF spectrum for the object, second part. Indicated hnes refer to Table 7. 
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Table 5. Identification of lines on the GBF spectrum. 



Identification 


A 


A fit 


Flux 


EW 


FWHM 




lim 




mJy 


A 


A 


Pa (8-3) 


0.9549 


0.9546 


9 


-21 


19 


Pa (7-3) 


1.0052 


1.0049 


10 


-15 


18 


He I (74)(3p3p0-6d3D) 


10314 


1.0315 


4 


-2 


15 


Fe II (z^Ft - b^Gg) 


1.0501 


1.0498 


34 


-24 


20 


?_ 


?_ 


1.0531 


6 


-6 


30 


He I (2s=^S - 2p=^P°) 


1.0832 


1.0830 


119 


-47 


15 


+He I (2s=^P - 2p3pO) 


1.0833 










?_ 


?_ 


1.0857 


40 


-30 


29 


He I (S^D - 6iF°) 


1.0917 


1.0911 


34 


-10 


12 


+ ?- 


?. 


1.0911 


34 


-10 


12 


Pa (6-3) 


1.0941 


1.0938 


110 


-59 


22 


Fe II (z^Fa - h^G^) 


1.1126 


1.1122 


46 


-17 


17 


?_ 


?_ 


1.1145 


7 


-3 


19 


He I (3p3po - 5d=^D) 


1.1970 


1.1970 


10 


-3 


27 


Nl (36)(3p2D3/2-3d2F7/2) 


1.2469 


1.2461 


25 


-11 


43 


He I (3s=^S - 4p3p0) 


1.2528 


1.2524 


19 


-3 


18 


[Fe II] (a^Dg - a^Dr) 


1.2567 


1.2568 


17 


-7 


47 


He I (3d3D - Sf^pO) 


1.2789 


1.2785 


49 


?-6 


14 


Pa (5-3) 


1.2822 


1.2817 


634 


-117 


22 


He I (3diD - 5fiF°) 


1.2976 


1.2979 


7 


?-l 


18 


I (3s3p - 48^8) 


1.3165 


1.3174 


-15 


2 


20 


?He I (3pip° - Ss^S) 


1.3415 


1.3427 


48 


-6 


21 


Br (25-4) 


1.4971 


1.4967 


21 


-1 


14 


Br (24-4) 


1.5005 


1.5003 


21 


-2 


21 


?Br (23-4) 


1.5043 


1.5034 


38 


-4 


26 


Br (22-4) 


1.5087 


1.5083 


49 


-5 


25 


+He I (4pipo - 3siS) 


1.5088 










Br (21-4) 


1.5137 


1.5135 


40 


-5 


29 


Br (20-4) 


1.5196 


1.5192 


36 


-4 


25 
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Table 5 — Continued 



Identification 


A 




Flux 


EW 


FWHM 








mJy 


A 


A 


Br (19-4) 


1.5265 


1.5259 


38 


-16 


21 


Br (18-4) 


1.5346 


1.5341 


65 


-6 


23 


Br (17-4) 


1.5443 


1.5439 


68 


-6 


22 


Br (16-4) 


1.5561 


1.5560 


87 


-12 


36 


Br (15-4) 


1.5705 


1.5699 


110 


-10 


25 


N III+C III 13-9 


1.575 


1.5753 


65 


-8 


33 


+Fe II (z2lii/2 - 3dHs%y2) 


1.576 










Br (14-4) 


1.5885 


1.5880 


147 


-14 


27 


[Fell] (a%/2 - a^Da/s) 


1.598 


1.5987 


5 


-1 


28 


Br (13-4) 


1.6114 


1.6107 


158 


-14 


27 


? 


? 


1.6168 


18 


-3 


51 


He II (25-8) 


1.6241 


1.6245 


12 


-1 


13 


Br (12-4) 


1.6412 


1.6407 


203 


-23 


37 
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Table 6. Identification of lines on the GRF spectrum, part 1. 



Identification 


A 




Flux 


EW 


FWHM 




//m 




mJy 


A 


A 


Br (18-4) 


1.5346 


1.5336 


39 


-4 


26 


Br (17-4) 


1.5443 


1.5432 


47 


-4 


25 


Br (16-4) 


1.5561 


1.5556 


66 


-10 


40 


Br (15-4) 


1.5705 


1.5694 


79 


-8 


27 


N III+C III 13-9 


1.575 


1.5747 


59 


-10 


48 


+Fe II (z2l?i/2 - 3d54s2ln/2) 


1.576 










Br (14-4) 


1.5885 


1.5877 


108 


-11 


30 


[Fe II] (a4F7/2 - a^Dg/s) 


1.598 


1.5950 


9 


7-6 


758 


Br (13-4) 


1.6114 


1.6104 


118 


-12 


32 


? 


? 


1.6159 


21 


-5 


70 


He II (25-8) 


1.6241 


1.6240 


9 


-37 


857 


Br (12-4) 


1.6412 


1.6408 


150 


-21 


44 


Br (11-4) 


1.6811 


1.6796 


217 


-33 


52 


Fe II (z^Fg/s - C^Fg/2) 

He I (4p3po - Ss^S) 


1.688 


1.6869 


328 


-32 


34 


1.7007 


1.7000 


118 


-9 


29 


[Fe II] (a^F5/2 - a^Dg/s) 


1.711 


1.7107 


?15 


7-1 


727 


Br (10-4) 


1.7367 


1.7371 


179 


-45 


90 


+ Fe II (z4D7/2 - C^P5/2) 












?_ 


?_ 


717406 


7160 


7-15 


748 


He II (19-8) 


1.7717 


1.7722 


717 


7-2 


732 


Br (9-4) 


1.8181 


1.8174 


473 


-36 


30 


Pa (4-3) 


1.8756 


1.8744 


597 


-22 


20 


He I (3pip° - 4diD) 


1.9094 


1.9093 


247 


-18 


31 


Br (8-4) 


1.9451 


1.9442 


386 


-29 


39 


Fe II (z^Fs - c^Fs) 


1.9746 


1.9736 


196 


-14 


36 


? 


? 


1.9856 


42 


-3 


35 


? 


? 


2.0026 


-156 


2 


8 
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Table 7. Identification of lines on tlie GRF spectrum, part 2. 



Identification 


A 




r lux 


w 


r WrIM 






//TTl 


m I^^ 


A 


A 


[Fc II] (a^ ,2 - a^Pg/s) 


2.046 


2.0457 


24 


-1 


26 


He I (2s^S - 2p^P°) 


2.0587 


2.0581 


628 


-42 


35 


Fe II (z^Fg/s - C^Fg/s) 

He I (3p=^P° - 4s=^S) 


2.089 


2.0881 


201 


-13 


36 


2.1126 


2.1116 


63 


-5 


43 


?N III + C III 


2.116 


72.1151 


7 42 


7 -4 


7 44 


?Mg II (5s2Si/2 - 5p2po 2) 


2.138 


2.1356 


135 


-8 


35 


?_ 


?_ 


2.1422 


71 


-6 


47 


+?Mg II (5s2Si/2 -5p2p?,2) 


2.144 










Br (7-4) 


2.1661 


2.1655 


698 


-45 


36 


He I (4pipo - 7diD) 


2.1847 


72.1857 


755 


7-7 


50 


Nal (4p2po „ - 4s2Si/2) 


2.2056,2.209 


2.2059 


45 


-3 


40 


[Felll fa^Hii/o - a^Go/o) 


2.224 


2.2240 


12 


-0 


21 


Fe II fz^D^/9 - c^P'^rt) 


2.240 


2.2393 


30 


-2 


40 


He I (4s^S - 6p^P°) 


2.3069 


2.3068 


17 


-1 


40 


He II (20-9) 


2.314 


2.3150 


20 


-1 


22 


?Pf (5-37) 


2.3218 


72.3215 


78 


7-0 


712 


Pf (5-36) 


2.3242 


2.3244 


15 


-1 


35 


+Fe II (z'^D^/o - C^Pq/o) 


2.3247 










-HC III f5D^P° - 5s^S') 


2.324 










'^Pf f5-32) 


2 3365 


2.3381 


25 


-1 


29 


-H'^Pf (5-31) 
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6. DISCUSSION 

6.1. Distance Estimate 

Simple qualitative arguments show that the source is likely to be galactic, and is not an 
outer galaxy : 

• the galactic latitude is —27' ; usually absorption prevents observations of extragalactic 
objects so close to the Galactic plane (however see Marti et al. (1998) and Ribo et 
al (2004) for two counterexamples) ; however, as INTEGRAL is aimed primarily within 
the galactic plane, a selection effect of unusual extragalactic objects is possible ; 

• the lines in the spectrum show no cosmological redshift ; relatively nearby objects (up 
to few Mpc) would show no redshift, but wc can exclude the possibility of a quasar ; on 
the other hand, the source is point-like, with a seeing better than 1", in disagreement 
with the possibility of a nearby AGN ; 

• in the X-rays, given the variability time scale, which leads to a limit of the distance 
between the X-ray source and the fluorescent emission region, the width of the Fe Ka 
fluorescent line is too small to come from a Seyfert II galaxy or an Ultra Luminous 
Infrared Galaxy (Revnivtsev et al. 2003; Walter et al. 2003). 

Estimating the distance from the r/D ratio of the SED fit of Figures 5 and 6 requires making 
some assumptions about the source. Wc assume that our source is a sgB[c] star, as we will 
discuss in subsection 6.3, and take the list of Lamers et al. (1998), giving the bolometric 
luminosity and the temperature for known sgB[e] stars. Using Stefan's law, we compute a 
star radius and use it to break the degeneracy of our fit. This gives a distance between 0.9 
(Hen S35 in LMC, type Bllab) and 6.2 kpc (R66 in LMC, type B8). As high mass stars 
remain usually close to star forming regions, we assume that our source is in a galactic arm. 
It can be then in the Sagittarius- Carina arm (0.7 kpc), the inner Scutum-Crux arm (3.2 kpc), 
in the Norma-Cygnus arm (4.8 kpc), and a star-forming complex at 7 kpc, where features 
and distances are taken from Russeil (2003). The next intercept between the direction of 
the source and a galactic arm, which is the Perseus arm at 10.8 kpc, is slightly too far. 
Conversely, all the stars of the list give a spectrum undistinguishable from the absorbed 
black body flt when put at the distance and absorption of our source, but it is only a simple 
consistency test of the method, not a clue that our source is really a sgB[e]. However, errors 
are difficult to assess, because of the various parameters involved, often in a non-linear way, 
and also because the characteristics of the sgB[e] stars show too much dispersion (with a 
possible existence of two subgroups, see Lamers et al. (1998)) to be considered as standard 
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Table 7— Continued 
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Fig. 10. — The observed integrated fluxes of hydrogen Uncs divided by the relative intensities 
from Hummer & Storey (1987). The lines indicate the best fit absorption model for the 
Brackett lines in the GBF spectrum (data in diamonds, fit in dots), in the GRF spectrum 
(triangles, dashes), all lines including Paschen (squares, solid line). Error bars are 4 % 
relative in flux, except for the lines at the edges and in the noisiest regions. 
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candles. 

Using the results of 4.2, and fixing the distance and the absorption, we can have an estimate 
of the luminosity. This is shown on Figure 11. For Ay = 17 A (hence T = 18 000 K), we 
have log(L/LQ) > 4 for D > 1 kpc. The corresponding radius is i?+ > IO-R0. These char- 
acteristics are compatible with a BO to B5 type giant star, or B or A supergiant star. For 
Ay = 11.8 for instance (hence T = 4 000K), we have log(L/L0) > 2 for D > Ikpc. The 
corresponding radius is i?^ > SO Rq. These characteristics are compatible with a G to M 
type of supergiant star. However we have seen that a higher absorption value is preferred 
by the data. 

A consistency check can be made by putting the source on a color-color HR diagram. 
This is done on Figure 12, using a HR diagram computed from template stars of Ruelas- 
Mayorga (1991), for the two absorption values considered above. The low absorption value 
is compatible with a distant red supergiant, or a close red giant. The high absorption case 
is compatible with an early-type supergiant, an OB or a B, provided that the distance is 
not greater than 4 kpc ; this suggests a location in the Scutum-Crux arm. The absorption 
cannot be substantially higher to be compatible with this figure ; this is consistent with 
Figure 6, showing that absorption is always below 18 for temperatures below 60 000 K. 

6.2. Nature of the Source 

The fact that IGR J16318— 4848 was detected during an outburst of X-ray emission, 
and that a previous outburst was found in the ASCA archival data of 1994 (Murakami et 
al. 2003), strongly points towards a binary system, including a compact object such as a neu- 
tron star or a black hole. A correlation between the X-ray and the radio emission in low/hard 
state black hole binaries is proposed in Gallo et al. (2003). According to this correlation, 
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Fig. 11. — Relation between distance and luminosity assuming Ay — 17 A (solid curve), and 
Av = 11.8 (dashed). 
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Fig. 12.— 
Mayorga 
and Ay — 



Position on the color-color HR diagram computed from template stars of Ruelas- 
(1991), for distance from 1 to 8 kpc, and for two absorption values : Ay — 17.4 
11.8. 
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and given the XMM flux and absorption in the 2-10 keV band (Matt & Guainazzi 2003), a 
radio flux of 43 ± 30 m Jy is expected (where the error bar comes from the uncertainties given 
in Gallo et al. (2003) for the parameters of their flt), more than one order of magnitude 
above the upper hmit given by ATCA (Walter et al. 2003). This is a hint that the compact 
object is a neutron star instead of a black hole, but given that the environment is so unique, 
one has to be cautious about applying this X-ray/radio relationship. Note that, similarly, 
the flux of CI Cam (12 mJy at 1.4 GHz according to Hjellming & Mioduszewski (1998)) is 
somewhat smaller that the expected flux (52 ± 30 mJy) given the ~ 65 mCrab flux during 
the outbursts (Prontera et al. 1998) if the correlation of Gallo et al. (2003) holds. 
We detect no signiflcant variation of the lines in wavelength, but given the short duration of 
the observations (10 minutes between the first and the last spectrum in each GBF and GRF 
bands), this is not surprising. Long duration observations would be useful in order to detect 
such a variation, if any, and to infer constraints on the masses of the system. 
Using the X-ray observations by XMM and INTEGRAL, Matt & Guainazzi (2003) and 
Walter et al. (2003) infer the following unabsorbed luminosities at 1 kpc : 1.3 x lO'^'^ crgs^^ 
between 2 and 10 keV, 3 x lO^^ergs"^ between 20 and 50 keV. As the photon index is very 
close to 2, the bolometric luminosity is somewhat difficult to estimate accurately, because it 
strongly model dependent, in particular of the cut-off energies. Assuming the X-ray emis- 
sion is dominating the SED of the compact object, we choose the sum of the two quoted 
values to be a good estimator of the luminosity (indeed it is a lower limit). Using the above 
estimation of the distance, this leads to a luminosity ranging from 1.3 x lO^^ergs"^ (0.9 
kpc) and 6.2 x lO'^^ergs^^ (6.2 kpc), at least 30 times lower than the Eddington limit of 
1.8 x lO^^ergs"^ for a neutron star of 1.4 M©. That suggests that moderate accretion is 
taking place, but given that the environment is so unique, one has to be cautious about 
applying this relationship. Note that the luminosity in the X-rays during the outburst is 
compatible with the common values for Be/X-ray binaries, where the compact object is a 
neutron star (Negueruela 1998). 

The huge difference in the absorption between the X-ray and our optical/NIR observations 
(the H I column density is 100 times higher in the X-rays than in our observations, see 
Walter et al. (2003)) indicates that the X-ray absorption is intrinsic to the compact ob- 
ject (Revnivtsev et al. 2003). It is probably due to mass loss from the companion star 
enshrouding the compact object, thereby absorbing the X-ray emission. Using the variabil- 
ity time scale of the X-ray emission, and the maximum delay observed between the Fe Ka 
fine and the continuum variations, Walter et al. (2003) deduce that the distance r between 
the zone in which fluorescent emission takes place and the X-ray source is constrained by 
r < 10^^ cm = 140 i?0. We suppose that the distance between the source of the NIR emis- 
sion and the source of the X-ray emission is at least r ; a hint for this hypothesis is that the 
absorption in the X-rays is much higher than in the NIR, although one has to be cautious. 
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because the nature of the absorbing material might not be the same in both wavelength 
domains. Unless the donor is a red giant of type K or later, this excludes a symbiotic binary 
system ; however, we have seen that our data prefer hot stars, and therefore that a High 
Mass X-Ray Binary (HMXB) is preferred. Assuming that we have a neutron star of 1.4 Mq 
orbiting a B type star with 1.4 < log(M/M0) < 1.7, and that we have r = 10^^ cm, we 
obtain a lower limit for the period between 28 and 38 days. With r = 10^^ cm, the minimum 
period is shortened to 0.9 to 1.2 days. In both cases, the period is far too long to have been 
detected by our observations. This point needs further observations to detect and measure 
an orbital period and would increase noticeably our knowledge of IGR J16318— 4848. It is 
possible that the X-ray outburst is close to the passage of the compact companion at the 
periastron, as it is the case in type I Be/X binaries (Negueruela 1998). 



6.3. About the Mass Donor 

The probably high absorption (we have seen that Ay — 17.4 is favored by our data) 
points to a massive supergiant star, with a luminosity higher than 10^"^ L© for a distance 
greater than 1 kpc, 10^ Lq for a distance greater than 2 kpc. The hypothesis of binarity leads 
to high mass X-ray binary (HMXB), the traditional mass donors of which being classical 
Be and OB supergiant stars. The latter, favored by the high luminosity, is not probable, 
because it generally leads to persistent X-ray emission. 

The fact that the lines of the spectrum are nearly all emission lines, and the fact that among 
them we have forbidden lines are key constraints for the identification of the spectral type 
of the mass donor. Usual stars have absorption hues. The P-Cygni profiles seem to be 
not numerous enough to make a classification as a P-Cygni star as HDE 316285 (HiUier et 
al. 1998). Bright infrared [Fe II] emission appears to be a common property among LBVs 
with prominent nebulae (Smith 2002) ; however, the line at 1.643 fim is very bright in these 
cases, whereas we do not observe this line, and the other forbidden lines are weak. Among 
the class of massive stars in transition, the spectral properties around 2 fim. are in good 
agreement with both B[e] and Of/WN stars (Morris ct al. 1996). The other star types (O, 
BO) considered in Morris et al. (1996) are disfavored because the Br7 line is in absorption 
in these types, but as noted in this article, there is always at least one example from each 
group that can be easily associated with another one, and therefore robust conclusions are 
difficult to infer from the infrared spectrum alone. However, Wolf-Rayet stars can be ruled 
out, as these stars are usually very poor in hydrogen. 

The NIR spectrum of IGR J16318— 4848 shows strong hydrogen Paschen, Brackett and Pfund 
emission lines. Although we have no spectroscopic data in the visible region, we can assume 
that Balmer emission lines are probably present at a strong level. The spectrum shows 
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also emission lines of Fe II, both allowed and forbidden transitions. The absorption seems 
to be higher than the estimated interstellar absorption along the line of sight (considering 
Ay = 11.8 ± 1.6 for interstellar absorption, or considering the very unusual colors of the 
source in respect to the field stars, as seen on Figure 4) ; moreover, the presence of Na I lines 
indicates the presence of circumstellar regions shielded from direct stellar radiation, but the 
absence of molecular H2 suggests a rather hot temperature. All this seems to converge to- 
wards the existence of hot circumstellar dust, although the lack of observations in the optical 
range makes difficult to assess an excess of IR emission. These points are the characteristics 
of the "B[e] phenomenon", according to the classification of Lamers et al. (1998). 
The fact that the luminosity is high suggests a supergiant B[e] (sgB[e]). Some of the sec- 
ondary criteria for a sgB[e] classification according to Lamers et al. (1998) are fulfilled : some 
lines show a P-Cygni profile indicating mass loss (criterion Bl of Lamers et al. (1998)), the 
possibly high He/H ratio indicates an evolved star (criterion B3), although this ratio is not 
very firmly constrained by our data, and a high extinction (criterion B4). However, we did 
not see any hybrid spectrum (there is no broad absorption feature) nor strong interstellar 
bands (contradicting criteria B2 and B4) , but the latter feature is only reported as " usual" . 
The non detection of the source in the USNO Bl.O catalog suggests that photometric vari- 
ations can be higher than 2 magnitudes in the R band (contradicting criterion B5, saying 
that the variation is on the order of 0.1 to 0.2 magnitude) ; however, the 2MASS magnitude 
in the J band (in 1999) is in agreement with our own measurements (2003). The criteria 
for the other B[e] subtypes (pre-main sequence stars, compact planetary nebulae or symbi- 
otic binaries) are clearly not satisfied. In addition, as already stated, the spectrum is very 
similar to that of CI Cam, both in the X-ray domain (Boirin et al. 2002), and in the NIR 
domain (Clark et al. 1999). CI Cam has been classified as a sgB[e] in the optical domain 
by Robinson, Ivans & Welsh (2002) and Hynes et al. (2002), as the first sgB[e] associated 
with a HMXB. The lack of IRAS detection is a concern with this hypothesis, as CI Cam is a 
significant IRAS source. However, CI Cam is also much brighter than the ICR J16318— 4848 
in the NIR and visible wavelength. By comparing the fiuxes in the K band of both sources, 
we expect a fiux at 12 ^m about 1 Jy or less for IGR J16318— 4848, around the upper limit we 
estimated for IRAS. Therefore, the lack of detection is not really conclusive. A classification 
as sgB[e] appears to be a reasonable choice, even if a classification as unclB[e]/sgB[e] may 
be more secure because of the problematic fulfillment of some secondary criteria. 
Therefore, we propose sgB[e] as a tentative identification for the massive companion of 
IGR J16318-4848, making of this source the second HMXB with a sgB[e] after CI Cam. 
There is evidence that the circumstellar material of a sgB[e] is concentrated in a disk (Hubert 
& Jaschek 1998), with a rarefied polar wind. As already proposed by Hynes et al. (2002) 
in the case of CI Cam, it is possible that the X-ray outburst is caused by the passage of the 
compact object in the disk, perhaps near the periastron. In this case, the burst should be 
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periodic, as the ASCA observations in 1994 suggest (Murakami et al. 2003). 



7. CONCLUSION 

The source IGR J16318-4848 was the first source discovered by INTEGRAL. In the 
course of a ToO program using the NTT telescope, we performed photometric and spectro- 
scopic observations less than one month after its discovery in the optical and NIR domains. 
We list here our main results : 

• we discovered the optical counterpart and confirmed an already proposed NIR candi- 
date ; 

• we performed an independent astrometry for this candidate ; 

• we obtained photometric measurements for the R, I and J bands, and flux upper limits 
for B and V, flux lower hmits for H and Kg ; 

• with the continua of our GBF and GRF spectra, these photometric measurements and 
with X-ray, radio and archival data, we constructed a SED covering 10 decades in 
wavelength ; 

• on the spectrum wc identified 72 emission lines, including forbidden lines, 80 % of these 
lines have been detected in CI Cam, suggesting a similar nature ; 

• the data favor an absorption of Ay — 17.4, greater than the interstellar absorption, 
and a temperature above 10 000 K ; 

• the data favor the existence of a dense circumstellar material, with stellar wind ; 

• the distance is between ~0.9 and 6.2 kpc ; 

• we propose as the most hkely hypothesis that the source is a HMXB with a sgB[e] star 
as the mass donor ; it would be the second case after CI Cam ; 

Complementary observations are needed in order to confirm our results, among them we 
propose : 

• high resolution NIR spectroscopy, if possible extended to optical, in order to : 

- extend the SED in the optical and in the mid-infrared to directly see a NIR excess ; 
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— check if the similarity with CI Cam as observed in the NIR is still valid in the 
optical ; 

— improve our results concerning P-Cygni profiles, and line broadening ; 
• long term follow-up spectroscopy and photometry, in order to : 

— seek for line variability ; 

— seek for a periodic behaviour to infer the orbital elements. 

This source shows many unusual features, the first is its strong intrinsic absorption. Inter- 
estingly, among the ten sources that INTEGRAL has discovered in this region, this feature 
is common (at least in the X-rays), to the three sources discussed by Revnivtsev (2003) : 
ICR J16318-4848, ICR J16320-4751 and ICR J16358-4726, although the A^h column den- 
sity is lower by an order of magnitude in the two latter systems (Rodriguez et al 2003; 
Patel et al 2004). However, a clear identification for the optical/NIR counterpart has been 
done only for ICR J16318— 4848. Moreover, the type of the mass donor, as inferred from our 
study, has been considered up to now as very rare. There is therefore the possibility that 
INTEGRAL, with the discovery of ICR J16318— 4848, has unveiled a new class of binaries 
that will deserve much attention in the future. 

This work uses observations made at the ESO NTT upon detection and locahzation 
made by INTEGRAL and XMM. We are very grateful to the ESO staff for their availability 
and skills for performing override programs, and special thanks go to Malvina Billeres for 
having performed these service ToO observations. We use the 2MASS, DENIS, USNO, 
GSC and DSS catalogs for astrometry and photometry. We acknowledge P. Ferrando, P. 
Goldoni, F. Lebrun, M. Ribo and J. Rodriguez for useful discussions and careful reading 
of the manuscript. SC thanks R. Hynes, I. Negueruela and M. Ribo for useful discussions 
on the nature of the source, and R. Walter for useful discussions. We thank the referee 
for a careful reading of the manuscript and useful suggestions. PF also acknowledge the 
CNRS/Federation de Recherche APC for funding. 
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